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Abstract 



Measurements of Z7* production are presented using data collected by the 
DELPHI detector at centre-of-mass energies ranging from 183 to 209 GeV, 
corresponding to an integrated luminosity of about 667 pb^^. The measure- 
ments cover a wide range of the possible final state four-fermion configurations: 
hadronic and leptonic {e'^e~qq, fi~^fi~qq, qqvv), fully leptonic {l^l^l^l^) and 
fully hadronic final states {qqqq, with a low mass qq pair). Measurements of 
the Z'j* cross-section for the various final states have been compared with the 
Standard Model expectations and found to be consistent within the errors. In 
addition, a total cross-section measurement of the l^l^l ^ cross-section is 
reported, and found to be in agreement with the prediction of the Standard 
Model. 
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1 Introduction 

The study of four-fermion processes in e^e~ interactions becomes increasingly im- 
portant as the centre-of-mass energy and the corresponding luminosity increase. The 
main goal of such studies is to verify the Standard Model predictions and to look for, 
or to set limits on, possible contributions arising from mechanisms beyond the Standard 
Model: for instance, anomalous triple gauge couplings [1] can usually give contributions 
to four-fermion final states. Moreover, such processes form an important background 
to new particle searches, such as those for charginos, neutralinos or non-standard Higgs 
bosons, and deviations from the Standard Model expectations would be a signal of new 
physics. LEP has provided a unique opportunity to study four-fermion interactions at 
several energies. On-shell pair production of W [2] and Z [3,4] bosons has been studied 
extensively. The focus of this paper is the measurement of the cross-section of neutral 
current processes with a Z and an off-shell photon (Z7* in the following). To this end, 
several channels were studied: l^l~qq (/ = e,/x), qquu, ^ {1,1 = e,ii,T) and qqqq 

(with a low mass qq pair). In addition, for I'^l'l +Z ~ final states, a measurement of the 
total neutral current cross-section has been performed. 

Figure 1 shows the main Feynman diagrams for four-fermion neutral current produc- 
tion in e~^e~ collisions. When there are no electrons in the final state these processes are 
dominated by the conversion processes shown in figure la. This graph represents eight 
different diagrams, usually referred to as the NCOS diagrams: two of them (usually re- 
ferred to as NC02) lead to ZZ production, two to 7*7* and four to Z7*. These four Z7* 
conversion diagrams are sketched in figure 2, and the square of their summed amplitudes 
is used in the definition of the signal to be measured in this paper, as explained below. 
A prominent feature of the graphs in figure 2 is the very different scale of the momentum 
transfer at the Zff and 7// vertices, an issue which must be properly addressed by the 
simulation programs (see section 3). For final states with electrons, other processes, such 
as t-channcl 7 exchange accompanied by Z*/7*-strahlung (figure Ic) and multipcriphcral 
production (figure Id), contribute significantly. In particular, the processes originating 
from the bremsstrahlung diagram (Ic), usually referred to as Zee and 7*ee, represent an 
important background to the measurements with electrons presented in this paper (see 
sections 4 and 6). Interference effects of these processes with those originating from Z7* 
can also be important and have to be taken into account. 

The Zj* production cross-section depends weakly on the centre-of-mass energy, but 
strongly on the mass of the virtual photon. For real photon production, e^e~ — > Z7, the 
cross-section reaches values above 100 pb, while in the kinematic region of Z7* production 
considered in this paper, its value is generally in the region of a fraction of a picobarn. 
Furthermore, in the Z7* production processes, particles coming from the conversion of 
low mass 7*s into hadrons or leptons are preferentially produced at very small angles 
with respect to the beam direction. A measurement of this cross-section has thus to be 
performed for a specific selection of the 7* mass and production polar angle. 

Data collected by the DELPHI experiment in 1997-2000 at centre-of-mass energies 
from 183 to 209 GeV were used, corresponding to an integrated luminosity of about 
667 pb~^. Results for each channel are given in the form of a comparison of the predicted 
numbers of selected events with those found in data. Combination of channel results into 
an overall Z7* cross-section is then performed. The resulting measurement is compared 
to the Standard Model expectation. The results presented here complement and augment 
those reported in previous studies of neutral current four-fermion production at LEP [5] . 
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1 a) Conversion 1 b) Annihilation 




Figure 1: The Feynman diagrams for four-fermion neutral current production in e^e 
collisions. 

This paper is organised as follows. Two definitions of the Z7* signal are given in 
section 2, one, the "Matrix Element definition", according to the Feynman diagrams 
contributing to Z7* production, the other, the "LEP definition", using invariant mass 
cuts. Short descriptions of the detector, of the available data sets and of the simulation 
programs used in the analyses are given in section 3. The subsequent sections provide 
descriptions of the analyses used for the first signal definition for each of the channels 
studied: l'^l~qq (section 4), qqvV (section 5), (section 6, where a total cross- 

section measurement is also presented), and qqqq (section 7). The results using the Matrix 
Element signal definition are presented in section 8, while the analyses and results using 
the LEP signal definition are described in section 9. Conclusions are given in section 10. 
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2 Signal definition 



Two different signal definitions were adopted in the analyses presented in this paper: 

• The Matrix Element definition: For each of the final states considered, the 
signal was first defined by applying the following kinematic selection on all charged 
fermions at generator level: 

|cos^j±| < 0.98 , 

where 9f± is the polar angle of the charged fermion with respect to the beam axis. 
Events with one or more charged fermions not fulfilling these selections were con- 
sidered as background. Then, for the surviving events, the signal was defined as the 
Z'-y* contribution coming from the four conversion diagrams shown in figure 2. This 
was achieved by weighting the events in the selected generator-level sample by the 
quantity 

\Mall? ' 

where M-z-y* and Mali are the matrix elements for Z7* and for all the graphs in 
figure 1, respectively. Analogously, using the same weighting technique, the compo- 



nents obtained by weighting events by the quantities ^"^f^ and 1 — 



- — 1^~^'^|2 were considered as background: these components represent, respectively, 
the contributions arising from non-Z7* four-fermion processes (including ZZ and 
7*7*, which are also produced via conversion diagrams) and from the interference 
effects between Z7* and non-Z7* graphs. Expected rates were thus computed using 
generated events weighted by the appropriate ratio. Efficiencies were defined from 
the simulated event samples as the ratio of selected weighted events over all weighted 
events. 

• The LEP definition: The second definition was agreed between the LEP Collabo- 
rations in order to combine results in a meaningful way. It is based on invariant mass 
cuts at generator level and explicitly avoids the difficult regions of low di-fermion 
masses. Depending on the final state, the following cuts were applied on invari- 
ant masses of fermion pairs and, where relevant, on lepton production polar angles: 
Mqq > 10 GeV/c^, Mi+i- > 5 GeV/c^, | cos^;±| < 0.95. Furthermore, it was required 
that only one fermion pair in the event had an invariant mass, Mf+f-, satisfying 
\Mf+f- — Mz\ < 2r^, where Mz and Tz are the nominal mass and width of the Z 
boson. Only the three dominant channels in the final result combination {iJi^ H'qq, 
e'^e~qq and qqi'P) were analysed using the LEP signal definition. 

In the rest of this paper, when not explicitly stated, it is implied that the Matrix Element 
signal definition is being used. 



3 Detector description and simulation 

A detailed description of the DELPHI detector and a review of its performance can be 
found in [6,7]. For LEP2 operations, the vertex detector was upgraded [8], and a set of 
scintillation counters was added to veto photons in blind regions of the electromagnetic 
calorimetry, at polar angles around 9 — 40°, 9 — 90° and 9 — 140°. 

The integrated luminosity of 666.7 pb~^ collected by the DELPHI detector at centre-of- 
mass energies from 182.7 to 209 GeV was used in the analysis. The luminosities collected 
at various centre-of-mass energies are shown in table 1. 
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Figure 2: Neutral current conversion diagrams for the Z7* process. 

During the year 2000, one sector (1/12) of the main tracking device, the Time Pro- 
jection Chamber (TPC), was inactive from the beginning of September to the end of 
data taking, which corresponded to about a quarter of the 2000 data sample. The effect 
was taken into account in the detector simulation and the corresponding small change of 
analysis sensitivity for this period was considered in the extraction of the cross-sections. 

Simulated events were produced with the DELPHI simulation program DELSIM [7] 
and then passed through the same reconstruction chain as the data. The generation of 
processes leading to four-fermion final states, mediated by charged and neutral currents, 
was done with WPHACT [9,10], interfaced to the PYTHIA [11] fragmentation and hadro- 
nisation model. For the charged current part, WPHACT incorporates the 0{a) Double 
Pole Approximation [12,13] radiative corrections to the doubly resonant WW production 
diagrams via a weighting technique, with the matrix elements provided by the YFSWW 
generator [14] . At a general level, WPHACT performs fully massive calculations all over 
the phase space, includes higher-order corrections and uses the package QEDPS [15] for 
initial state radiation. Two additional features, particularly relevant for the analyses 
described in this paper, were implemented in WPHACT: the study of the most suitable 
scale to use for aqEo at the 7* vertices of the diagrams in figure 2, and the treatment of 
the hadronisation of low mass virtual photons. The first of these problems was solved in 
WPHACT by implementing the running of aqEo at the level of the event generation, thus 
using the value of the coupling constant corresponding to the mass of the photon propaga- 
tor at the 7* vertices. The second problem was addressed by interfacing WPHACT with 
a special package [16] for the specific treatment of the hadronisation of low mass qq sys- 
tems. This package provides a description of the hadronisation from the 7* — > qq process 
in the mass region below 2 GeV/c^ both due to the presence of hadronic resonances (with 
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Year 


1— 


Integrated 




[GeV] 


luminosity [pb~^] 


1997 


182.7 


55.0 


1998 


188.6 


158.1 


iyyy 


191.0 


zo.U 


1999 


195.5 


77.0 


1999 


199.5 


82.0 


1999 


201.6 


41.0 


2000 


205.0 


81.3 


2000 


206.5 


147.3 


Total 


197.1 


666.7 



Table 1: Luminosity- weighted centre- of- mass energies and integrated luminosities of the 
data analysed. During the year 2000, the energies reached were in the range 202-209 GeV 
and clustered mainly around 205 and 206.5 GeV. 

subsequent decays described by PYTHIA) and in the continuum, based on experimental 
e+e" data at low energy. This is particularly important for the qquv (section 5) and qqqq 
(section 7) analyses, which explore the low mass qq region. Phase space cuts are applied 
in WPHACT and can be found in table 2 of [10]. The study of the backgrounds due to 
qq{'j), ji^jji^ij) and T+r~(7) production was made using the KK2f [17] model; e^e^{^'i) 
events were simulated with BHWIDE [18]. Two-photon interactions were generated us- 
ing WPHACT for the regions in which the multiperipheral contribution is not dominant 
and using BDK [19] for the pure two-photon region; PYTHIA 6.143 was used to model 
two-photon processes with single and doubly resolved photons. 

4 Study of the l^l~qq final state 

In this section the analysis of the final state containing jets and a pair of identified 
muons or electrons is described. The two final state leptons in the process e^e~ l~^l~qq 
are usually well isolated from all other particles. This property can be used to select such 
events with high efficiency in both the muon and electron channels. Events with t'^t~ 
pairs have not been considered here. This part of the analysis follows very closely the 
one performed in [3], where an identical final state was studied. 

A loose hadronic preselection was first applied, requiring that the events have at 
least 7 charged particles and a charged energy above 0.30y^. To suppress the radiative 
return to the Z (final state on-shell Z production with the emission of a hard initial 
state radiation (ISR) photon) the event was rejected if a photon with energy more than 
60 GeV was found. The selection procedures were then carried out in a similar way for the 
fji^ fJTqq and e^e~qq channels. In order to maximise the lepton identification efficiency, 
any charged particle with momentum exceeding 5 GeV/c was considered as a possible 
lepton candidate around which nearby photons, if present, could be clustered. This 
was found to be necessary to improve the energy evaluation in the presence of final state 
radiation from electrons. In the case of the e^e^qq channel, photons with energy between 
20 GeV and 60 GeV were also considered as electron candidates in order to recover events 
in which the electron track was not reconstructed. For both electrons and muons, "strong" 
and "soft" identification criteria were then defined. Muons were considered as strongly 
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identified if selected by tlie standard DELPHI muon identification package [7], based 
mainly on finding associated hits in the muon chambers. For soft muon identification, 
only kinematic and calorimetric criteria were used. Electrons were considered as strongly 
identified when the energy deposited in the electromagnetic calorimeter exceeded 60% of 
the cluster energy or 15 GeV, whichever was greater, and when the energy deposited in 
the hadron calorimeter was less than a specified limit. For soft electron identification, only 
requirements on the momentum of the charged particle in the cluster and on the energy 
deposited in the hadron calorimeter were imposed. Electron candidates originating from 
applying the clustering procedure around a photon were considered as softly identified. 
Events with at least two lepton candidates of the same flavour and opposite charge were 
then selected.^ All particles except the lepton candidates were clustered into two jets and 
a kinematic fit requiring four-momentum conservation was applied, after appropriately 
adjusting the errors on lepton energies in cases where photons had been added by the 
clustering procedure. 

At least one of the two lepton candidates was required to satisfy strong lepton identi- 
fication criteria, while softer requirements were specified for the second. 

Two discriminating variables were then used for final event selection: p™*". the lesser 
of the transverse momenta of the lepton candidates with respect to the nearest jet, and 
the per degree of freedom of the kinematic fit. Cuts on these variables were applied, 
with values depending on the final state and on the quality of the lepton identification 
(see [3]). 

4.1 Results 

The distribution of the reconstructed mass of one fermion pair when the mass of the 
second one is within 15 GeV/c^ of the nominal Z mass is compared with the predictions of 
the Standard Model in figure 3, separately for j.!^ ji'qq and e^e~qq events. In the fJi^ ^~qq 
channel the Z7* contribution is clearly separated from the background component and 
is mostly concentrated in the region of the decay Z — > qq, as expected. In the e^e~qq 
case, the distribution is fiatter, indicating the presence of non-resonant diagrams. In both 
channels there is good overall agreement between the observed and predicted numbers 
of events. In the e^e^qq channel an accumulation of events is observed in the invariant 
mass distribution of the e^e~ pair in the region between 50 and 60 GeV/c^, with 7 events 
observed where 2.4 are expected. Various studies and comparisons with results of the 
other LEP experiments were performed, leading to the conclusion that this excess is most 
probably due to a statistical fluctuation. 

The bidimensional distributions in the plane of the masses of the two fermion pairs 
predicted by the Standard Model are shown in figure 4 for the two channels studied, sepa- 
rately for Z7* and background. The presence of non-resonant contributions, particularly 
of the type Zee and 7*ee, is clearly visible in the e'^e~qq case. The distributions were 
binned as shown graphically in figure 4, using a small number of irregularly sized bins. 
This allowed the regions where most of the background is concentrated to be avoided, 
except for the Zee contribution, while keeping as much signal as possible. Bin sizes were 
chosen in order to have an approximately equiprobable distribution of simulated events, 
with a finer binning in e~^e~qq so as to follow better the more complicated structure of 
the background distribution. The observed and predicted numbers of events selected by 
this procedure at each energy point are reported in table 2. 

^The requirement of having leptons of opposite charge was dropped in the case of candidate electrons originally identified 
as photons, for which no charge information is available. 
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Figure 3: l~^l~qq final state: Distributions of tlie mass of one fermion pair wlien tlie mass 
of the second is within 15 GeV/c^ of Mz- The two lower plots are for the e~^e~qq channel 
and the two upper plots for the fi~^fJ'~qq channel. The points are the data summed over all 
energy points, the dark (red) histogram is the distribution of the background predicted 
by the Standard Model, and the light (light blue) histogram is the predicted distribution 
of the Z7* signal. 

On combining the data from all energy points, the efficiency and purity of the selected 
/i+yU~gg sample were estimated from the simulation to be 42.0% and 84.7%, respectively. 
The background is composed of fi'^fi~qq events outside the signal definition and of contri- 
butions from other final states. Interference effects in the fJ'~^fi~qq channel in the region 
considered are negligible, as they account for less than 0.1% of the Z7* cross-section. 
The predicted composition of the background is shown in table 3. 

In the e~^e~qq channel, the purity of the selected sample is estimated to be only 49.2%, 
mostly because of the unavoidable Zee background, while the efficiency was evaluated 
to be 24.3%. Interference effects between Z7* and other four-fermion processes were 
estimated to account for -15% and are thus not negligible: they are mostly concentrated 
in the region of Z7* - Zee overlap. The predicted composition of the background is shown 
in table 3. 
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Figure 4: l^l~qq final state: Bidimensional distributions in the plane of the di-fermion 
masses predicted by the Standard Model for signal Z7* (left-hand plots) and background 
(right-hand plots) for the two channels studied, averaged over all energy points. The bins 
used for the fit are also shown. 



In order to disentangle the Z7* from the Zee contribution more effectively, the distri- 
bution of the polar angle of the direction of the e+e~ pair was studied as a function of the 
reconstructed invariant mass Mee in the range defined by the first five bins in figure 4. 
Correlation plots are shown in figure 5 for signal and background: the distributions are 
well separated because in the Zee case, which dominates the background, even e'^e~ pairs 
of large invariant mass are emitted at low polar angles, due to the t-channel nature of the 
production process. The binning in figure 4 was therefore modified, as shown in figure 5, 
by doubling each bin, depending on whether a) the polar angle of the direction of the 
e+e~ pair was in the barrel region (40° < 9ee < 140°) or in the endcap region {9ee < 40° 
or 6ee > 140°) for bins 1-5, and b) the polar angle of the direction of the hadronic system 
was in the barrel or in the endcap region for bins 6-7. A total of 14 bins was thus used 
for the e^e^qq cross-section measurement. This procedure resulted in an 8% reduction 
of the statistical error compared to the case where only mass bins were used. 

The value of the Z7* cross-section at each energy point was extracted using a binned 
likehhood fit technique (see section 8) and the values were then combined to get global 
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3.6 


3.8 


Total 


35 


36.7 


31.1 
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39 
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Table 2: Observed numbers of events in the qq and e^e qq channels at each energy 
compared with the Standard Model predictions for signal and background. 



Background source 


11+ fi qq 


e+e qq 


WW 


0.8 


1.6 


qq{i) 


0.1 


1.8 


T^T~qq 


2.5 


2.6 


non-Z7* l^l~qq 


2.2 


18.0 


Interference 


< 0.001 


-5.6 


Total 


5.6 


18.4 



Table 3: Composition of the background to Z^j* production in the fi+fi^qq and e+e^qq 
final states predicted by the Standard Model. The entries show the expected numbers 
of events, summed over all energy points. The row labelled non-Z7* l'^l~qq shows the 
four-fermion neutral current contributions from processes leading to the same final state 
as the signal, but defined as background, as described in section 1. 

results, separately for ii+^~qq and e'^e'qq. Only the value of the Z7* contribution was 
varied in the fit, while all non-Z7* contributions, backgrounds and interference terms 
were fixed to the Standard Model expectations. Figure 6 compares the data in each bin 
used in the fits to the (J^+n^qq and e+e^qq final states with the results of the fit, showing 
the contributions from the Z7* signal, from the non-Z7* component of each of the final 
states, from the interference terms, and from the other sources of background. 

These results were used to derive the combined values of the Z7* cross-section for the 
Matrix Element signal definition, as described in section 8. 

In the e^e~qq case, where the presence of non- resonant diagrams is relevant, a two- 
parameter fit was also performed as a consistency check, leaving both the Z7* and the 
non-Z7* contributions free to vary, while fixing the remaining background sources and 
interference terms to the Standard Model expectations. No significant change in the 
Z7* cross-section result was observed, while the ratio Rnon-z-y* of the measured to the 
predicted cross-section of the non-Z7* contribution to e'^e~qq was determined to be 
Rnon-z-y* = 1-15 lo!23) whcrc the error is statistical only. 
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Figure 5: Distributions of the mass of the electron-positron pair in e~^e~qq events for val- 
ues of Mee less than 80 GeV/c^ versus the polar direction of the pair for signal Zj* (upper 
plot) and background (lower plot) when the mass of the hadronic system is compatible 
with the Z mass. The plot shows the predictions of the Standard Model, averaged over 
all energy points. The binning adopted for these events follows that in figure 4 with an 
additional division into barrel and endcap regions, described in detail in the text. 



4.2 Systematic errors 

Several sources of systematic error were investigated. 

Uncertainties in lepton identification, signal efficiency and background levels were 
evaluated using a procedure similar to that in [3], where the same final states were 
studied. 

Uncertainties in the lepton identification were estimated by comparing semileptonic 
WW events selected in data and simulation using the strong lepton identification criteria. 
Uncertainties in signal efficiencies were evaluated by comparing the P™" and distribu- 
tions in data and simulation for all llqq events selected without mass cuts. Corresponding 
uncertainties in background levels were evaluated by comparing samples of events selected 
in data and in simulation, requiring both isolated tracks not to be identified as leptons, 
while maintaining all the other criteria. Finally, uncertainties in the background level in 
the e~^e~qq channel from fake electrons were studied with qq{'~f) events selected in data 
and in simulation with purely kinematic criteria. These effects and the statistical uncer- 
tainty of simulated data yielded a combined relative systematic error on the efficiency to 
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Figure 6: Results for fi'^ft'qq (upper plot) and e^e~qq (lower plots). The points are 
data, summed over all energy points; the shaded histograms show the results of fits to 
the Zj* component and the predicted Standard Model contributions from other sources. 
See figures 4 and 5 and the main text for the bin definitions. 
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select /^"'"/i qq and e+e qq events of ±5.0%, and a relative uncertainty in the background 
level of ±15%.2 

Systematic effects coming from the fitting procedure were investigated. Fit results 
were found to be stable within the expected statistical uncertainties against variations 
of bin sizes, the number of bins and, for e~^e~qq, the number of fitted parameters. No 
systematic error was thus attributed to this source. 

Possible systematic effects arising from the treatment in the fit of the Zj* interfer- 
ence term with the other contributions (particularly the non-Z7* one) were taken into 
account (see section 4.1) for e^e~qq. Both the one-parameter and two-parameter fits 
were repeated and the interference term was weighted with a factor proportional to the 
product of the Z^* and non-Z7* amphtudes. This changed the cross-section result by 
2%. Note, however, that this procedure neglects a possible change in the phase between 
the two interfering amplitudes with respect to that predicted in the Standard Model, and 
the procedure adopted therefore estimates the maximum possible effect that the unknown 
phase could have. A systematic uncertainty of ±2% was thus ascribed from this source 
for e'^e~qq events. 

The systematic error coming from the uncertainty in the luminosity measurement was 
evaluated to be ±0.6% both for e^e~qq and for fi~^fi~qq. 

The total estimated systematic errors on the measured Zj* cross-sections were ±5% 
for fi^fi~qq and ±6% for e'^e~qq. 

5 Study of the qqui? final state 

The qquP channel is observed in a final state topology of hadronic matter and sub- 
stantial missing energy. About half of the Z'-f* cross-section in this channel comes from 
the region of qq masses below 6 GeV/c^. Thus, final states often have the characteristic 
signature of "monojets" , with the low invariant-mass hadronic system, which is the event 
visible mass, arising from the 7* hadronisation and recoiling against a highly energetic 

pair which escapes detection. 

Three analyses were performed and combined. The first analysis was intended to probe 
the low mass region of the hadronic system, so as to be efficient in the region of virtual 
photon mass. My*, below 6 GeV/c^, where most of the cross-section is expected. It is 
denoted as the "low mass analysis" in the following. The second analysis exploited the 
large energy imbalance of qqiyp events, and retained some efficiency in the very low mass 
region of the hadronic system. It is denoted as the "energy asymmetry analysis" in the 
following. The third analysis was intended to have good overall efficiency for high My. 
at the expense of having very small efficiency in the low M^* region. It is denoted as the 
"high mass analysis" in the following. 

A common event preselection was defined for the three analyses, aimed mainly at 
reducing the backgrounds from two-photon and Bhabha events. The energy measured in 
the electromagnetic calorimeters was required to be less than 60 GeV in total and less 
than 10 GeV at polar angles below 15° and above 165°. Events with identified electrons 
at polar angles below 15° and above 165° were excluded; the visible energy of the event 
was required to exceed 15% of the centre-of-mass energy; the polar angle of the direction 
of the event missing momentum was required to be in the range 15° < Omiss < 165°; and 
at least two charged particles with momentum greater than 200 MeV/c were required. 

*^In both cases determinations were limited in a<;curax;y by the statistics of the available samples, and should be inter- 
preted as upper bounds. 
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An extensive use of veto counters was implemented in all three analyses: events with 
hits in the photon veto counters far from energy deposits in calorimeters or reconstructed 
tracks were rejected. The details of the algorithms adopted are given in the following 
sections. 

In order to increase the available statistics, no explicit lower cut on the reconstructed 
mass of the hadronic system was applied. 

The numerical values of the cuts applied to kinematic variables in the three analyses 
were chosen using an optimisation procedure described in section 5.4 below. 

5.1 Low mass analysis 

Events with a visible mass M^ig < 6 GeV/c^ and with visible energy E^ig larger 
than 20% of the centre-of-mass energy were selected. In addition, in order to limit the 
background from leptonic decays of l^s {W ^ e/fi u, W ^ rz/, r e/^ u), it was 
required that no identified muon be present, while at most one electron was allowed in 
the event and its energy was required to be less than 30 GeV. Furthermore, events with 
the polar angle of the direction of the missing momentum in the range 38° to 42° (which 
is insufficiently covered by calorimeters, see section 3) were rejected. The event was 
then split into two hemispheres by the plane perpendicular to the thrust axis: events 
were rejected if there were hits in the photon veto counters in the hemisphere containing 
the direction of the missing momentum, while events with hits in the veto counters in 
the opposite hemisphere were accepted only if their angular separation from the closest 
charged-particle track or calorimetric energy deposit was less than 20°. 

When used alone, this analysis selected 10 events in data and 6.7 in the simulation, of 
which 4.3 were signal and 2.4 were background. 

5.2 Energy asymmetry analysis 

In this analysis events were required to show a marked imbalance in the spatial distri- 
bution of the detected reaction products. Only events with total visible energy exceeding 
20% of the centre-of-mass energy were accepted. Then two hemispheres were defined by 
a plane perpendicular to the direction of the thrust axis, and the total energy in each 
hemisphere was estimated from the curvature of charged-particle tracks and from calori- 
metric measurements. It was required that the energy in one of the two hemispheres 
account for at least 99% of the total energy in the event. This was the main topological 
selection of the analysis and provided an implicit upper cut-off on the total visible mass 
of events. 

Signals from photon veto counters were used to discard events with possible loss of 
energetic photons in the insensitive regions of the electromagnetic calorimetry by adopting 

the same algorithm as in the low mass analysis (sec section 5.1). In order to limit further 
the background from processes which have most of the cross-section in the forward region 
(mainly Bhabha and two-photon events), the cut on the polar angle of the direction of 
the missing momentum was tightened and required to lie in the range 25° < 6 < 155°. 

At this level, the background was completely dominated by the WW and Wei/ pro- 
cesses. In order to reject leptonic decays of Ws, events with identified muons were 
discarded, while events with at most one electron were kept if the energy of the electron 
did not exceed 25 GeV and the electron was not isolated, i.e. its angle with respect to 
the closest charged-particle track was not larger than 10°. 
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Additional selections were implemented in order to suppress further the WW and Wei^ 
backgrounds. Part of this background arises from hadronic decays of one W, accompanied 
by undetected leptons coming from the decay of the other W or lost in the beam pipe 
(especially in the case of Weu). Such events usually show larger visible mass than is 
expected from signal events, due to the sizeable mass of the W boson. A selection on the 
event visible mass was thus imposed, requiring M^^g < 45 GeV/c^. Another important 
fraction of the remaining background comes from WW events with both l^s decaying to 
Ts, W Tfr-i with the visible decay products boosted into the same hemisphere. The 
signature of these events is that a few particles carry most of the visible energy and have 
visible mass above a few GeV/c^. Two more selections were imposed to reject such a 
source of background. Events with visible mass above 5 GeV/c^ and with more than 
90% of the visible energy carried by the two most energetic particles were discarded. The 
remaining events were forced into two jets with the LUCLUS algorithm [20]. Events with 
total particle multiplicity below 11 and an angle between the two jets above 30° were 
rejected. 

When used alone, this analysis selected 25 events in data and 29.5 in the simulation, of 
which 17.3 were signal and 12.2 were background. Half of the background was contributed 
by the WW and Wev processes. 



5.3 High mass analysis 

In this analysis a cut on the multiplicity of chargcd-particle tracks was applied, re- 
quiring it to be larger than 4. This implied that the efficiency of the analysis dropped 
essentially to zero for qq masses below 2 GeV/c^. The main topological selections were 
applied at jet level. Jets were reconstructed using the LUCLUS algorithm and the events 
were forced into a two-jet configuration. An upper cut on the opening angle of the two 
jets was set at 78°. The parameter dH^^^ was defined to be the value for which the event 
passes from a two-jet to a single jet configuration: only events with ^2"" < 30 GeV/c 
were retained. The acoplanarity (defined as the complement of the angle between the 
jets projected on the plane perpendicular to the beams) was required to be larger than 
90°. Then the event was split into two hemispheres about a plane perpendicular to the 
thrust axis and the energy asymmetry, evaluated as in section 5.2, was required to be 
larger than 95%. Events with missing mass less than 80 GeV/c^ were rejected. 

A further selection was imposed on the energy of the visible system, E^isi rejecting 
events with E^s > 80 GeV. In the absence of initial- and final-state radiation, the energy 
and the mass of the qq system in the Zj* process are related in the following way: 



2v^ 



The quantity Ekin = ^ ^2^"'' defined, using the visible mass of the event. It was 
then required that the difference between E^i^ ^-nd the visible energy of the event E^^g 
did not exceed 45 GcV. This cut, and the cut on Eyig described above, were effective in 
suppressing the WW and qq{'~f) backgrounds. 

Events with hits in the photon veto counters were accepted if the angular distance 
between these hits and the direction of the closest jet was less than 30°; otherwise they 
were rejected. 

When used alone, this analysis selected 21 events in data and 20.7 in the simulation, 
of which 13.4 were signal and 7.3 were background. Most of the background is due to 
WW and Weu events. 
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5.4 Results 

The three analyses were combined on an event-by-event basis, by selecting events 
which passed any of the three selections. Numerical values of the cuts were optimised 
in a two-stage procedure. First, for each analysis separately, all the cuts relevant to 
that analysis were varied such that the product of efficiency and purity of the selected 
sample was maximised. Then the most important cuts in each analysis were allowed 
to vary simultaneously, keeping the other cuts at the values obtained in the first stage, 
and the product of the efficiency and purity of the sample selected by any of the three 
analyses was maximised. (The values listed in sections 5.1, 5.2 and 5.3 are the result 
of this last optimisation procedure). In total, 42 events were found in data and 41.3 
in the simulation (with a total overlap between the three selections of about 30%); of 
the simulated sample, 23.4 events were signal and 17.9 were background. The most 
abundant source of background was predicted to come from Weu events, which accounted 
for 7.9 events, mainly in the channel qqeu and partially in Ti/ev. On-shell WW processes 
contributed about 4 events to the background, with 2.9 of them containing at least one W 
decaying to ru. The remaining main sources of background were qq (about 2 events), rr 
(about 2 events) and other four-fermion neutral current processes (1.1 events). Table 4 
shows the numbers of signal and background events predicted by the Standard Model 
and the observed numbers of events in the qquP channel at the various centre-of-mass 
energies. 



E (GeV) 


Data 


Total MC 


Signal 


Background 


182.7 


3 


3.5 


2.3 


1.2 


188.6 
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10.1 


6.0 


4.1 


191.6 
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1.3 


0.9 


0.4 


195.5 


7 


4.4 


2.9 


1.5 


199.5 


6 


5.2 


2.9 


2.3 


201.6 


2 


2.5 


1.3 


1.2 


205.0 


9 


4.9 


2.6 


2.3 


206.5 


5 


9.4 


4.5 


4.9 


Total 


42 


41.3 


23.4 


17.9 



Table 4: Observed numbers of events in the qquu channel at each energy compared with 
the Standard Model predictions for signal and background. 

The differential efficiencies of the three analyses as a function of the generated mass, 
M{qq), estimated from the simulation, are shown in figure 7, together with the efficiency 
for the combined selection. The overall selection efficiency, averaged over all masses, was 
estimated to be 38.8%. The distribution of the reconstructed visible mass, M^j^, for the 
42 data events is shown in figure 8, which also shows the distributions for the simulated 
signal and background events. Good agreement is observed with the Standard Model 
expectations. 

The value of the Zj* cross-section at each energy point was extracted using a count- 
ing technique and the values were then combined to get a global result. All non-Z7* 
contributions, backgrounds and interference terms were fixed to the Standard Model ex- 
pectations. The result was used to derive a combined value for the Z^* cross-section in 
the Matrix Element signal definition, as described in section 8. 
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Figure 7: Selection efficiency of the qqvv analyses, averaged over all energy points, as 
a function of the generated M{qq) mass. The efficiency is shown for each of the three 
analyses (see text) separately, and for the combined analysis. 

5.5 Systematic errors 

Various sources of systematic error were considered. 

The predicted background cross-sections were varied according to the following values: 
WW : ±2%, qq : ±5%, Weu : ±5%, rr : ±5%, four-fermion neutral current processes: 
±5%. The combined effect on the cross-section measurement was estimated to amount to 
±2%, with the main contribution coming from the uncertainty on the Weu cross-section. 

Uncertainties on the signal efficiency coming from the Monte Carlo generator were 
studied by comparing different generator models. In particular, a sample of qquu was 
generated with the EXCALIBUR [21] four-fermion generator for masses of the hadronic 
system M{qq) > 10 GeV/c^. For generated masses below 10 GeV/c^ the hadronisation 
model in EXCALIBUR is not as reliable as that in WPHACT and systematic effects 
from that region were evaluated separately (see below). The full analysis was applied to 
the EXCALIBUR sample and a difference of 3% in the signal efficiency was obtained. A 
systematic uncertainty of ±3% was thus conservatively ascribed to this source. 

Systematic uncertainties due to the description of the hadronisation mechanism in the 
qq system were taken into account. It was assumed that these effects can be relevant 
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Figure 8: Distribution of the visible invariant mass of the hadronic system in the qqvv 
selection, compared with the Standard Model predictions for signal and background. 
The points are the data, summed over all energy points, the light (light blue) histogram 
shows the predicted signal contribution, and the dark (red) histogram shows the predicted 
background. 

for masses M{qq) < 10 GeV/c^ (see above), affecting the analysis mainly through corre- 
sponding uncertainties in charged-particle multiplicity distributions. These effects were 
not expected to be large because two of the three analyses (the low mass analysis and the 
energy asymmetry analysis) adopted a very low cut on the charged-particle multiplicity. 
The study of these effects was split into two parts, corresponding to the resonance and 
the continuum contributions respectively (see section 3). In the simulated sample the 
dominant resonances were identified, their corresponding detection efficiencies computed, 
and their contributions varied by amounts derived from the uncertainties in their known 
measured cross-sections: 1% for p and production, 10% for resonances decaying to final 
states with 3 or 4 charged particles, and 30% for resonances decaying to states with 5 or 
6 charged particles. The effect on the estimated cross-section was found to be negligible; 
this is not surprising as p production, for which the cross-section is accurately deter- 
mined, accounts for about 80% of the cross-section below 2 GeV/c^. As a second step, 
the contribution of the resonances was subtracted from the hadronic mass distribution. 
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the charged-particle multiplicity distribution of the remaining sample studied and the 
analysis efficiency evaluated as a function of the number of charged particles. The effect 
of a possible error in the determination of the charged-particle multiplicity distribution 
was then estimated by stretching the observed distribution by +20%, rebinning, and 
applying the efficiency curve to the new distribution. The procedure was repeated, com- 
pressing the distribution by 20%. The range of cross-sections obtained from the stretched 
and compressed distributions was taken as an estimate of the systematic error. The effect 
on the Z7* cross-section amounted to ±4%. 

Another source of systematic uncertainty considered was the reliability of the simula- 
tion in correctly estimating the amount of background. As explained in section 5.4, the 
main backgrounds are WW events, with one or both l^s decaying to r, and Weu events, 
with the on-shell W decaying hadronically or to tu. These events share the common 
feature of having the decay products of one W detected on one side, and missing energy 
on the other side. The missing energy is due to the low angle electron, typically lost in the 
beam pipe in the Weu case, or to an undetected decay lepton or charged-particle track in 
the WW case. Furthermore, in both topologies, additional missing energy is carried by 
the escaping neutrino. In order to evaluate the reliability of the simulation in estimating 
the efficiency to detect backgrounds in such a topology, events with features similar to 
those of the background in the qquu analysis were studied. In particular, WW events 
with one W decaying to a detected lepton (electron or muon) or to an isolated charged- 
particle track, which was then artificially removed from the event, can mimic most of the 
WW and Weu background, with the second W playing the role of the hadronic signal. 
Therefore events with an isolated electron, muon or other charged-particle track were se- 
lected. Identified leptons or other charged-particle tracks were initially required to have 
momentum larger than 10 GeV/c and an angle with respect to the closest charged-particle 
track larger than 10°. The selected candidate track was then excluded from the event 
and the selections in sections 5.1, 5.2 and 5.3 applied to the remaining system. At the 
end of the procedure, 142 events with an isolated muon were found in data and 135.7 in 
the simulation, 110 events with an isolated electron were found in data and 115.1 in the 
simulation, and 79 events with a single isolated charged-particle track were found in data 
and 72.4 in the simulation. The distribution of the isolation angle of the selected lepton 
or single track after all the cuts is shown in figure 9. Good agreement between data 
and simulation is observed. The dominant contributions to the events selected in this 
way come from semileptonic WW production, rr events and, to a lesser extent, Bhabha 
events, and can thus be used to emulate the background to the qqi'i' signal: when the 
isolated lepton or other charged-particle track is excluded from the sample, the remaining 
system is strongly asymmetric in the angular distribution of the visible momentum, and 
of the same topology as the background expected in the qqvP sample. (The estimated 
contribution from the qqi/P signal to this sample is totally negligible). 

The agreement between data and simulation in figure 9 was subjected to more detailed 
checks, for example by selecting the region of the distribution in the isolation angle of 
the single charged-particle track which enriches the sample in WW events: two-fermion 
events preferentially populate the region of large isolation angle, being almost back to 
back, and their contribution can be greatly reduced with a cut at around 130°. Similarly, 
other checks were made for different visible-mass and track-multiplicity regions; in all 
the cases the agreement between data and simulation was good within the errors. The 
statistical error in the total of 331 events selected by this procedure was thus taken as 
an estimate of the systematic uncertainty due to the background evaluation from the 
simulation; this gave a contribution of ±3% on the qquP cross-section measurement. 
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Figure 9: Isolation angle of selected electrons, muons and single charged-particle tracks 
in the sample selected to mimic the background to the selected qqvv sample. The points 
are the data, summed over all energy points, the light (green) histogram is the predicted 
WW contribution, and the dark (red) histogram is the rest of the background. 

Systematic uncertainties from the trigger efficiency were investigated and found to 
be negligible: the triggering efficiency for a single charged-particle track with transverse 
momentum > 3 GeV/c is already very well determined [22], while in the present 
analyses a charged-particle track multiplicity of at least 2 was required, with transverse 
momenta of selected events in general well in excess of 3 GeV/c. 

The systematic uncertainty coming from the luminosity measurement was estimated 
to give an error of ±0.6% on the cross-section measurement. 

The statistical error from the limited simulated sample gave an uncertainty of ±5%. 

Finally, the stability of the result as a function of the applied experimental cuts was 
checked by varying the selections of the three analyses, first separately and then at the 
same time. The procedure set up to maximise the product of the efficiency and purity of 
the simulated sample (see section 5.4) was used to vary all the relevant cuts of sections 5.1, 
5.2 and 5.3 within reasonable values; selections were accepted if the predicted value for the 
selected sample differed by less than the statistical error of the optimum value obtained 
from the simulated sample used in the analysis. For each such selection, the background 
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level and number of events in data were estimated and a value for the qqiyp cross-section 
was measured. The root mean square of the distribution of the cross-sections evaluated 
in this way was estimated to amount to 3%. As this number is compatible with the 
statistical fluctuations intrinsic to this procedure, no systematic error was added. 

The total estimated systematic error on the Z^* qquP cross-section measurement was 
thus estimated to be ±8%. 

6 Study of the l^l~t^t~ final state 

The Feynman diagrams of figure 1 give rise to six possible final states with four charged 

leptons: fi'^ fj,~ fj,~^ fi~ , e~^e~e~^e~, t~^t~t~^t~ , e~^e~fi'^fi~, e~^e~T~^T~ and fi'^fi~T~^T~ . 
These final states have a rather clean experimental signature, but do not contribute 
significantly to the total four-fermion production cross-section due to the low branching 
fraction of Z/Y ■ 

The selection of events in the ~ final state was restricted to topologies with 

four well reconstructed charged particles with momenta larger than 2 GcV/c (henceforth 
called lepton candidates). Events with two additional well-measured charged particles 
with opposite charges were allowed, provided that the pair was compatible with a pho- 
ton conversion, or that the momentum of both particles was less than 2 GeV/c. Five 
additional charged particles were allowed in the event if their tracks did not point to the 
vertex; such tracks were not considered in the following steps of the analysis. The previous 
selections implied that for e+e^r^r~, /i^/i.~r+r~and t^t~t^t~ events only one-prong r 
decays were considered. The sum of the charges of the lepton candidates had to equal 
zero and the angle between the directions of any two of them had to be larger than 5°. 

The four lepton candidates were required to fulfil the following additional selection 
criteria: the momenta of at least three of them had to exceed 6 GeV/c, their total energy 
had to be greater than 0.25\/s (to reject background from two- photon interactions), and 
the length of at least three of the candidates' tracks was required to be greater than 50 cm. 
Beam-gas and r+r~7 events were rejected by requiring that the four lepton candidates 
were not all in the same hemisphere with respect to the beam direction. For data taken 
during 2000, in the period when one sector of the TPC was not working, a slightly more 
relaxed criterion for track selection was applied if the track traversed that sector. 

Selected events in the data were compared with simulated signal and background 
samples generated at the eight centre- of- mass energies. The expected numbers of events 
for signal and background, together with the numbers of events found in data, are shown 
in table 5 both for the full sample of I'^l't'^l'" events, selected as described above, and 
for the ^7* sample defined in section 6.2. The overall +Z ~ selection efficiency is 
~15%, increasing slightly with ^/s for the full sample, while for the Z7* selection it 
ranges between 22% and 30%. The most important contribution to the non-/"*"/"/ ~ 
background comes from e^e" — > e'^e~qq events with low qq mass. The second most 
important contribution is due to the e'^e" — > T~^T~{'y) process. Good agreement was 
found between the data and the predictions of the simulation after each selection was 
applied sequentially. For the Z7* sample, the main background is due to ~ 
contributions from non-Z7* processes. 
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t 1 7 7' J_ 7' 

z+/ / +/ 


" full sample 




t 1 7 1' 1 -1^ 


Z7* sample 




Data 


Total MC 


Signal 


Background 


Data 


Total MC 


Signal 


Background 


182.7 


3 


3.9 


3.4 


0.4 


1 


1.5 


0.5 


1.0 


188.6 


14 


12.4 


10.0 


A 

2A 


2 


4.8 


1.6 


3.2 




i 


1 
1.0 


1.6 


0.2 


i 


O.D 


0.2 


0.4 


195.5 


5 


5.3 


4.6 
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2.0 


0.7 


1.3 
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2.7 


2.4 


0.4 


2 


0.9 


0.3 


0.6 


205.0 


7 


5.3 


4.8 


0.5 


4 


1.8 


0.7 


1.2 


206.5 


7 


9.6 


8.2 


1.4 


3 


3.1 


0.9 


2.1 


Total 


48 


47.0 


40.1 


6.8 


17 


16.6 


5.5 


11.1 



Table 5: Observed numbers of events in the channel for the full sample and for 

the Z7* sample at each energy, compared with the Standard Model predictions for signal 
and background. In the case of the Z7* sample, the background contributions arc defined 
to include the non-Z7* '^l ~ contribution and the non-l~^l~l '^l ~ contribution. 



6.1 Particle identification and final state classification 

Events selected in the final state were classified into one of the six final 

states according to the number of identified muons, electrons and pions. A constrained 

fit procedure was also used to complete the identification. 

Muon identification was performed by combining the standard DELPHI identifica- 
tion package [7] in the muon chambers with the energy deposition profile in the hadron 
calorimeter and the energy deposited in the electromagnetic calorimeter. 

Electron identification required that there be no signal in the muon chambers and 
no energy deposited in the hadron calorimeter after the first layer. The energy in the 
electromagnetic calorimeter in a 2° cone surrounding the candidate particle was required 
to be larger than 1 GeV. For electrons satisfying these criteria, the momentum of the 
charged particle was replaced by the energy deposited in the electromagnetic calorimeter. 

Pions were identified as tracks leaving an energy deposit in the electromagnetic 
calorimeter compatible with a minimum ionizing signal, no hits in the muon chambers 
and energy deposited in the layers of the hadron calorimeter compatible with the profile 
of a hadron shower. 

The assignment of the final state proceeded as follows: 

• If no e'^e" or /i^/x" pair was identified, the four particles were considered as r decays 
and the final state to be t~^t~t~^t~; 

• If two pairs were identified as e+e", or r+r", the final state was considered 
to be fully identified; 

• If one e'^e~ or pair was identified and the second pair had two identified 
particles, different from one another, the event was considered to be e'*"e~T"^T~ or 
ji^ li~ T~ , respectively. The second pair was also designated as t^t~ if only one 
particle was identified and was different from the identified pair, or if neither was 
identified; 

• If the event had 3 identified electrons or muons and one unidentified particle, two hy- 
potheses were considered: that the 4 particles were identical or that the unidentified 
particle was one of a r+r" pair. 
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A constrained kinematic fit was then performed on tlie selected events, imposing four- 
momentum conservation. This imphes a four-constraint fit in the case where both lepton 
pairs are either electrons or muons, and a two- or zero-constraint fit in the cases where, 
respectively, one or two tan pairs are assumed present, as the magnitude of the tau 
momentum was taken to be unknown. Where more than one kinematic hypothesis could 
be applied to the same event, the decision procedure and the final identification were 
based on the probabihty of the of the fit and the relative errors of the fitted masses. 
In the case of four identical particles, the combination for which a pair of leptons had 
reconstructed mass within 15 GeV/c^ of the nominal Z mass was chosen or, if this 
condition was not fulfilled, the combination with the largest fitted invariant mass of 
a pair of leptons was selected. If no acceptable hypothesis was found, further fits were 
tried where kinematically possible, assuming, in addition to the four leptons, the presence 
of an unobserved ISR photon in the beam pipe; again the best resulting fit was selected. 
Figure 10 shows the distributions of the larger and smaller mass pairs for the full data 
sample, calculated for each event from the results of the chosen fit, and compares them 
with the predictions of the Standard Model. 

The efficiencies for assigning the correct final states to the selected events were es- 
timated from the simulation. The results are summarised in table 6, which shows the 
expected numbers of events from the full "'"/ ~ sample which were identified in each 
of the possible final states, as well as the efficiency and purity. 



Identified 
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T+T T+T 


Total 


Bck 


e+e~e+e~ 
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0.0 
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0.0 
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92 
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0.0 
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0.0 


0.4 


1.6 


0.0 


6.0 


0.4 


7 
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0.4 


0.5 


0.6 


0.0 


0.2 


0.2 


1.9 


1.5 


2 


6 


Efficiency (%) 


10 


11 


11 


13 


12 


9 





Table 6: Upper six rows: Expected numbers of signal and background events and purity 
for each identified final state for the full ~ event sample, estimated from the 

simulation. The number of events found in the experimental data is also given for each 
final state. Bottom row: Estimated efficiency for selection and correct classification of 
each state with respect to the total content of the sample. 

In 2% of the cases the events could not be classified in any of the six final states, as there 
was no complete a priori identification of all the particles in the event and the constrained 
fit failed. Due to lack of identification of electrons or muons, mainly in regions with poor 
coverage by the electromagnetic calorimetry or muon chambers, or from inefficiencies 
in the particle identification algorithms, a substantial fraction of '^l ~ events was 
misidentified as having a pair of taus. The 48 events selected in the data were classified as 
follows: 7 in the e+e^e+e^ channel, 14 as e'^ e" /i'^ /i" , 16 as e+e~T"'"T~, 2 as /i'^ /i" /i'^ /i" , 
7 as ii'^ii~T~^T~ and 2 as t~^t~t~^t~ . 

6.2 ZY production in Results 

The value of the cross-section at each energy point was extracted using a procedure 
which followed closely that adopted for the l^l~qq channels in section 4.1. Bidimensional 
mass distributions were constructed in the plane of the masses of the pairs with the larger 
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Figure 10: Four-lepton channel: Fitted masses of two selected lepton pairs: bigger mass 
(top plot), smaller mass (bottom plot), compared with Standard Model predictions. The 
points are the data, summed over all energy points, and the histograms represent the 
predicted contributions to the selected event sample. In the legend, "background" means 
the contribution from non-/"*"/"/ "*"/ ~ final states. 
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and smaller mass in the event. The distributions were binned using the same definition 
as the first five bins in figure 4 for the e~^e~qq case. Of the 15.1 events predicted as 
the Z7* plus the non-Z7* contributions to the total signal, 3.2 were predicted 

in the e^e~e^e~ channel, 4.5 in e~^e~n'^fi~, 3.0 in e~^e~T'^T^ . 1.6 in fi^ fjr fj,^ fjr . 2.2 in 
fi~^fi~T^T~ and 0.6 in t~^t~t^t~ , while, of the 17 selected data events, 1, 4, 6, 1, 5 and 
were assigned to each of these channels, respectively. 

A one-parameter binned likehhood fit to the Z7* contribution was per- 

formed, fixing the non-Zj* contribution and the remaining backgrounds and interference 
terms to the Standard Model expectations. These results, shown in the right side of 
table 5, were used to derive the combined values of the Zj* cross-section in the Matrix 
Element signal definition, as described in section 8. 

6.3 Z7* production in l^l~t^t~'. Systematic errors 

Several sources of systematic uncertainties were investigated. 

The main contribution to the systematic error in the track selection came from the 
difference between data and simulation in the number of reconstructed charged-particle 
tracks. In order to estimate this uncertainty, samples of dimuon events were generated 
and the numbers of events with one, two or three reconstructed charged tracks com- 
pared in data and simulation. Prom the comparison, a conservative uncertainty of ±5% 
was assigned as the systematic error from this source. For dimuon events with two re- 
constructed charged-particle tracks, the difference between data and simulation in the 
number of events with total charge equal to zero was found to be of the order of +0.5%. 

A contribution of 1.5% was added due to differences between data and simulation in 
the charge misidentification of electrons in the low polar angle region. 

Systematic uncertainties originating from particle identification were also taken into 
account. Two pure samples of e^e~ and fi^fi^ final states were selected from the data 
using particle identification criteria independent of those described in section 6.1 and were 
compared with simulated samples of the same final states. Then the identification criteria 
for electrons and muons were applied to both samples and the difference in the efficiencies 
between data and simulation was taken as a systematic error. This resulted in errors of 
±0.5% for muons and ±5% for electrons. The poorer of the two estimates was also used 
for taus and adopted as a systematic uncertainty on the cross-section measurement. 

Possible errors arising from the procedure adopted in the fits to the mass distri- 
bution were studied. Several checks were performed, in close analogy to those described 
in section 4.2. First, simulated samples of events with electrons in the final state (which 
receive large contributions from t-channel processes) were split into two categories, de- 
pending on whether or not the electrons were identified in the event reconstruction. 
The cross-sections of the two samples were measured and then combined. Secondly, a 
one-parameter fit to the mass distribution was performed, both on the whole selected 
^ sample and on the two separated samples with final state electrons described 
above, allowing only the Z7* component to vary. From the spread of the results of these 
additional fits, a systematic error of ±7% was estimated. 

The error in the efficiency for selecting signal events due to the limited Monte Carlo 
statistics was evaluated to be ±0.6%. The limited statistics available for the different 
background processes were also taken into account, as well as the theoretical uncertainties 
in the cross-sections, resulting in contributions of ±0.06% and ±1.1%, respectively. Fi- 
nally, a contribution to the systematic error of ±0.6% was estimated from the uncertainty 
in the measurement of the luminosity. 



25 



The total estimated systematic error on the Z7* cross-section measurement 

was thus estimated to be ±10%. 

6.4 Measurement of the total cross-section for produc- 
tion 

In this section we report a total cross-section measurement for production, 
in addition to the study of ZY production in the four-lepton topology described in 
section 6.2 above. 

As the cross-section does not vary too much within the energy range of LEP2, all the 
data and the Monte Carlo simulations for the different energies were grouped together. 
The total cross-section was then estimated from a likelihood fit to the Poissonian proba- 
bility for observing the number of events found in the data, given the expected number 
corresponding to a total cross-section, a, for l^l~l +Z ~ production, plus the estimated 
number of background events (see table 5). 

The total cross-section for the ~ processes was found to be 

a = (0.430 ± 0.072 ± 0.023) pb 

within the visible region, defined by | cos^^l < 0.98, at a luminosity- averaged centre-of- 
mass energy of 197.1 GeV. The first error quoted is statistical; the second is the estimated 
systematic error, derived as described in section 6.3 above, but without including effects 
involving particle identification. 

This result is in good agreement with the predicted cross-sections from WPHACT, 
which range from 0.440 pb at ^/s — 182.7 GeV to 0.375 pb at a/s = 206.5 GeV, giving a 
luminosity-weighted average cross-section of 0.403 pb within the visible region at = 
197.1 GeV. 

7 Study of the qqqq final state 

The measurement of the Z7* contribution in the qqqq channel has to deal with back- 
ground processes such as 95(7) and WW which have cross-sections larger by orders of 
magnitude than the signal. It is thus not feasible to measure the Z7* cross-section in all 
the possible qq mass spectrum. Only a restricted region was thus considered here, for 
low values of the reconstructed mass of one qq pair. The signature of the process studied 
in this analysis is the presence of a highly energetic isolated low mass jet from the 7* 
hadronisation (preferentially directed in the forward region), recoiling against a system 
of two (or more) jets from the hadronic Z decay. The study of the 7* system was limited 
to final states with only two charged particles and an arbitrary number of neutral parti- 
cles; this choice was driven by the expectation that, in the low mass region, the process 
7* — » qq is dominated by the hadronisation chain 7* ^ p° vr+vr". Furthermore, an 
explicit cut on the reconstructed mass of the two selected chargcd-particle tracks was 
used, as explained below. The Z7* signal definition was kept the same as in the other 
channels studied (with no limits on the 7* mass); as a result, the two selection criteria 
mentioned above (those requiring low charged-particle multiplicity and low reconstructed 
mass) imply a large inefficiency in the analysis of events with 7* — > qq for 7* masses above 
a few GcV/c2. 

The principal backgrounds arise from production of ^^(7) , WW and final states from 
other four-fermion neutral current processes such as qq/j,'^/!'', qqe^e" and qqT'^T~ . 
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A pre-selection was applied to the data in order to select hadronic events compatible 
with the expected topologies. The total charged-particle multiplicity was required to be 
larger than 20; the ratio \^ had to be larger than 77%, where is the reconstructed 
effective centre- of- mass energy [23]; events with neutral particles with electromagnetic 
energy exceeding 50 GeV were excluded; the missing energy of the event was required 
to be less than 82% of the centre-of-mass energy; and the number of identified muons 
was required to be less than two (to limit the background from qqn^ ti~ events). Events 
were then clustered according to the LUCLUS [20] algorithm with the parameter djoin 
set to 6.5 GeV/c, and it was required that the number of reconstructed jets in the event 
be larger than two. One of the jets had to contain at least one charged particle with 
momentum exceeding 32 GeV/c and to have charged-particle multiplicity of two, while 
an arbitrary number of neutral particles was accepted in the jet. The pair of charged 
particles was then subjected to the selections listed below: 

• The impact parameters of the two charged particles were required to be compatible 
with production at the primary event vertex; 

• The total energy of the pair was required to be larger than 63 GeV; 

• The two charged particles had to be of opposite charge; 

• The total energy deposited by the two particles in the electromagnetic calorimeters 
was required to be less than 40% of the total energy of the pair; 

• Identified muons and electrons (soft identification criteria, see section 4) were not 
allowed in the pair; 

• The system recoihng against the jet containing the selected pair was forced into 
a two-jet configuration and the full kinematics of the three jets was completely 
determined by their space directions. Then the two-jet system not containing the 
selected pair was required to have a reconstructed mass within 11 GeV/c^ of the 
nominal Z mass; 

• The invariant mass of the two charged particles had to be less than 2.1 GeV/c^. 

Numerical values of the cuts were optimised by scanning the full range of the relevant 
discriminating variables and calculating, for each set of values, the cross-section and the 
product of the efficiency and purity of the selected sample. The set with the highest 
value of the product of efficiency, e, and purity, p, corresponding to e = 2.2% and p = 
69.6%, was chosen, yielding a ratio signal ^ _ rpj^^ procedure selected 7 events 

a/ background 

in data and 6.9 in the simulation, of which 4.8 were signal and 2.1 were background. 
The main backgrounds came from WW (1.1 events), ^'^(7) (0.4 events) and other four- 
fermion neutral current processes (0.4 events). Figure 11 compares the distribution of 
the reconstructed mass of the pair of selected charged-particle tracks before the last cut 
with Standard Model predictions. Table 7 shows the predicted numbers of signal and 
background events and the observed numbers of events in the qqqq channel at the various 
centre-of-mass energies. 

The value of the Z7* cross-section at each energy point was extracted using a counting 
technique and the vahics were then combined to determine a global result. All non-Z'y* 
contributions, backgrounds and interference terms were fixed to the Standard Model 
expectations. The results were used to derive a combined value for the Zj* cross-section 
for the Matrix Element signal definition, as described in section 8. 

7.1 Systematic errors 

Various sources of systematic error were considered. 
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Table 7: Observed numbers of events in the qqqq channel at each energy compared with 
the Standard Model predictions for signal and background. 

The predicted background contributions from WW, qq{j) and four-fcrmion neutral 
current production were varied by changing the cross-sections for these processes accord- 
ing to the values given in section 5.5: the combined effect on the cross-section measure- 
ment was estimated to amount to ±0.8%. 

The statistical error corresponding to the limited simulated sample gave an uncertainty 
of ±8%. 

The reliability of the simulation in reproducing the amount of background was checked 
by repeating the analysis, selecting pairs of particles of the same charge. The same cuts 
as those described in section 7 were applied, with the exclusion of the requirement on the 
total charge of the pair. No events were selected in data, while 0.56 were predicted by the 
simulation. The results are of course compatible, but to derive a numerical estimate for 
a systematic error, the procedure was modified so as to select a larger number of events: 
the cut on the invariant mass of the pair of chargcd-particle tracks - made at 2.1 GeV/c^ 
in the main analysis - was increased to 10 GeV/c^. All the other selections were left 
unchanged. This gave 3 events in data and 4.5 in the simulation, of which 3.6 were due 
to WW production and 0.5 to qq{'y) backgrounds. 

A similar study was performed to check the four-fermion neutral current background, 
which gave a negligible contribution in the previous procedure. The selections in section 7 
were repeated on data and simulation, but replacing the veto on identified electrons or 
muons in the selected pair of charged tracks by the requirement that at least one of the 
two tracks was positively identified as a lepton (electron or muon). In addition, the cut on 
the invariant mass of the pair was softened to 10 GeV/c^, as for the check described in the 
previous paragraph. This resulted in 8 events selected in the data and 6.9 predicted from 
the simulation, of which 5.7 were due to the four-fermion neutral current background (in 
particular I'^l'qq events, with I = e,jjL, r) and 0.9 from the WW background. 

As the two last procedures (requirement on the total charge of the pair and on the 
presence of leptons in the pair) each showed good agreement between data and the pre- 
dictions of the simulation, the results were summed, and the larger of the statistical error 
of the data and the difference between data and simulation was assumed as a systematic 
uncertainty. This was estimated to be ±13% on the cross-section measurement. 

The uncertainty on the cross-section measurement due to the measurement of the 
luminosity was evaluated to be ±0.6%. 

Finally, the stability of the result as a function of the applied experimental cuts was 
checked by varying the numerical values of the analysis selections. The procedure set 
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Figure 11: The distribution of the reconstructed invariant mass of the selected pair of 
charged-particle tracks in the qqqq analysis, compared with the predictions of the Stan- 
dard Model. The points are the data, summed over all energy points and shown before 
the final selection of Mpair < 2.1 GeV/c^; the light (blue) histogram shows the predicted 
Z7* contribution, and the dark (red) histogram shows the predicted background. 



up to maximise the product of the efficiency and purity of the simulated sample (see 
section 7) was used to vary all the relevant cuts within reasonable limits; selections were 
accepted if the predicted product of the efficiency and purity of the sample differed by 
less than the statistical error of the simulated sample from the optimum value used in the 
analysis. For each new selection, the signal efficiency, background level and number of 
events in data were estimated, and a value for the cross-section was measured. The root 
mean square of the distribution of the cross-sections thus obtained was evaluated to be 
15% of the central value. As this number is compatible with the statistical fluctuations 
intrinsic to this procedure, no systematic error was added. 

The total estimated systematic error on the qqqq Z^* cross-section measurement was 
thus estimated to be ±15%. 
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8 Results 

The measurements described in the previous sections all show good agreement with 
the expectations of the Standard Model. In this section, we use these measurements 
to give results for the ratio, Rz^y*, of the measured to the expected Z'-f* cross-section, 
for each of the final states considered, for their combination at each of the LEP energy 
points at which data were taken, and for the overall average. All these results are given 
in terms of the Matrix Element signal definition (see section 2). Results for the LEP 
signal definition are given in section 9. 

Individual cross-sections were extracted by maximising probability functions with re- 
spect to the value of the Zj* cross-section: Poissonian probabilities, based on the number 
of events selected in data and predicted in the simulation, were used for the qquP and qqqq 
channels; probability functions derived from fitting procedures were used for the fJ>'^fJ^~qq, 
e^e~qq and /"'"/"/ "channels. For each centre-of-mass energy, results were expressed in 
terms of the ratio Rz-y* of measured to expected cross-sections, thus automatically taking 
into account the (smooth) dependence with energy predicted by the Standard Model. 
The results obtained for the different energies were first combined for each channel sepa- 
rately, and then into a single value. Global likelihoods were constructed to perform such 
combinations. The central value was defined as the point of minimum — log L distribu- 
tion and the statistical error as the interval around the central value which contained 
68.27% of the probabihty. The results obtained for the different channels are shown in 
table 8 and in figure 12. The table also shows the average cross-section predicted by the 
Standard Model for each of the final states considered at the luminosity-weighted average 
centre-of-mass energy of 197.1 GeV. 

Table 9 compares the results at the various energy points, averaged over the different 
channels, with the Standard Model predictions, and this comparison is also shown in 
figure 13. 

The systematic uncertainties for each channel were studied by introducing appropri- 
ately modified assumptions for backgrounds and efficiencies (as described in the corre- 
sponding sections) and were considered as fully correlated between the energies. The 

effect of systematic uncertainties in the combination of different channels was taken into 
account considering the uncertainties due to the luminosity measurement and to varia- 
tions in the predicted background cross-sections as correlated between the channels, and 
all other effects as uncorrelated. 
The final result is 

i?Z7* = 1.04 lo:i2(stat) ± 0.04(syst) 

for I cos6'/±| < 0.98, as shown in tables 8 and 9 and in figure 12. This result is in good 
agreement with the Standard Model expectation. 

9 Analyses and Results for the LEP signal definition 



The analyses of the three dominant channels in the final result combination {^^fJ^^qq, 
e^e~qq, qqvP, see table 8), described in sections 4 and 5, were repeated adopting the 
LEP signal definition (see section 2). Some modifications were introduced to the analyses 
described in the sections referred to above, in order to take into account the fact that 
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channel 


Rz-y* 


a(pb) 


^SM (pb) 




0.98+^:^^ ± 0.05 


0.1081^:^^^ ±0.005 


0.11 


e^e~qq 


i.05i^::^g ± 0.06 


0.115l^:^i± 0.007 


0.11 


qqvv 


1.051^:^^ ± 0.08 


0.084^^:^1? ± 0.006 


0.08 


l+l-l'+l'- 


i.3i;^i^±o.i3 


0.039l^:^J|^ ± 0.004 


0.03 


qqqq 


1.09;n7±0.16 


0.316^^:1;^^ ± 0.047 


0.29 


Total 


1.04l^:{^ ± 0.04 


0.6661^:^?^ ± 0.026 


0.64 



Table 8: Ratios of measured to predicted cross-sections and measured cross-sections for 
individual channels contributing to the Z'y* process, using the Matrix Element signal 
definition (see section 2). The first errors are statistical and the second systematic. In 
the last column ^^^/(pb) is the average, luminosity- weighted Z7* cross-section predicted 
by the Standard Model at the average energy of 197.1 GeV. 



E (GeV) 


Rzj* 


aiph) 


(pb) 


182.7 


1.551^;::;^ ±0.01 


i.i5i;!:^^ ± 0.03 


0.71 


188.6 


0.831^:^^ ± 0.04 


0.57t'^:'^i ± 0.03 


0.69 


191.6 


0.41^^:?? ± 0.04 


0.27^^:^^ ± 0.03 


0.67 


195.5 


1.18^^:1^ ± 0.04 


0.78^^:^^ ± 0.03 


0.66 


199.5 


0.891^:1^ ± 0.04 


0.58l^:i ± 0.03 


0.65 


201.6 


2.631^:?^ ± 0.04 


1.661^:^^ ± 0.03 


0.63 


205.0 


1.52^^:^^ ± 0.04 


0.90^^:i ± 0.02 


0.59 


206.5 


0.44t^:^^ ± 0.04 


0.25l^:it ± 0.02 


0.57 


Average 


1.04l^:i^ ± 0.04 


0.671^:^^ ± 0.03 


0.64 



Table 9: Ratios of measured to predicted cross-sections and measured cross-sections av- 
eraged over the different channels at the various energy points, using the Matrix Element 
signal definition (see section 2). The first errors are statistical and the second systematic. 
The last column shows the Standard Model predictions, and the entries in the last row 
refer to the luminosity- averaged centre-of-mass energy of 197.1 GeV. 

the di-fermion invariant mass regions below the cuts described in section 2 must now be 
considered as background. 

• In the l^l^qq analysis two additional selections were introduced with respect to those 
described in section 4: it was required that the reconstructed mass, M;+;-, of the 
two charged leptons be larger than 4 GeV/c^ and that the reconstructed mass of the 
remaining hadronic system be larger than 8 GeV/c^. This corresponds to reducing 
the content of bin 1 in the plots of figures 4, 5 and 6 and of bin 4 for muons and 
bin 6 for electrons in the same figures. The other steps of the analysis were left 
unchanged and the same procedures were applied to evaluate the systematic errors. 
The total systematic uncertainty on the measured Z7* cross-section with the LEP 
signal definition was estimated to be ±6% for lJi~qq and ±7% for e~^e~qq. 

• In the qquu analysis similar modifications were introduced. The low mass analysis 
(see section 5.1) was not used, while in the energy asymmetry and high mass analyses 
(see sections 5.2 and 5.3, respectively), it was required that the reconstructed mass 
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Measured/Predicted Z7* Cross— Section 

DELPHI 



eeqq 

qqqq 



o.gslfg ± 0.05 

1 .o5:;jo ± 0.06 

1 .05!;i? ± 0.08 

^.5^T4l± 0.13 

1.09!o°:4'?± 0.16 



Average 1 .04I12 ± 0.04 



Figure 12: Ratios of measured to predicted cross-sections for individual channels con- 
tributing to the Z'j* process, using the Matrix Element signal definition (see section 2). 
The vertical band displays the total error on the combination of the channels. 



of the hadronic system be larger than 8 GeV/c^. The other steps of the analyses 
were left unchanged and the same procedures were applied to evaluate the systematic 
errors. The total systematic uncertainty on the measured Zj* cross-section with the 
LEP signal definition was estimated to be ±16%. 

The same procedures as described in section 8 were applied in order to obtain results 
for the Z'J* cross-sections with the LEP signal definition. The final results for the three 
channels used are summarised in table 10. A combined value of 

azr = 0.136 to:oi(stat) ± 0.008{syst) pb 

was obtained for the luminosity-weighted cross-section with the LEP signal definition, in 
good agreement with the Standard Model prediction of 0.151 pb. 
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Figure 13: Combined Zj* cross-section as a function of the centre-of-mass energy, using 
the Matrix Element signal definition (see section 2). The solid line is the Standard Model 
prediction; the dashed line represents a 5% uncertainty around this prediction. The full 
(red) point is the average cross-section result, plotted at the luminosity- weighted average 
centre-of-mass energy. 



10 Conclusions 

In the data sample collected by the DELPHI detector at centre-of-mass energies rang- 
ing from 183 GeV to 209 GeV, the values of the Zj* cross-section contributing to the 
four-fermion final states fi^fi~qq, e~^e~qq, qquu, l^l^l^l^ and qqqq with | cos6'j±| < 0.98 
have been measured and compared with Standard Model expectations. A combined value 
of 

Rz^* = 1.04 toA2istat) ± OM{syst) 

was obtained for the ratio of the measured to the predicted cross-section in the Matrix 
Element signal definition (described in section 2). This corresponds to a luminosity- 
weighted measured cross-section of 

azr = 0.666 tomristat) ± 0.026{syst) pb , 

in good agreement with the value of 0.640 pb predicted by the Standard Model. 
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channel 




a(pb) 


^SM (pb) 




0.74^^:^^ ± 0.05 


0.031^^:^t^ ± 0.002 


0.042 


e^e~qq 


1.051^:^^ ±0.08 


0.061^^:^{{ ± 0.004 


0.058 


qqvv 


0.83i^;^|±0.13 


0.042^^:^^^ ± 0.007 


0.051 


Total 


o.90l^:i^ ± 0.05 


0.1361^:^^? ± 0.008 


0.151 



Table 10: Ratios of measured to predicted cross-sections and luminosity-weighted cross- 
sections for individual channels contributing to the Z7* process, using the LEP signal 
definition (see section 2). The first errors are statistical and the second systematic. In 
the last column cj^m (pb) is the average, luminosity- weighted Z7* cross-section predicted 
by the Standard Model. 

Additional cross-section measurements in the channels n'^n~qq, e^e'qq and qquv were 
performed using the common LEP signal definition (also described in section 2). A 
combined, luminosity- weighted, value of 

az-y" = 0.136 tQ^{stat) ± 0.008(s?/st) pb 

was obtained, in good agreement with the Standard Model prediction of 0.151 pb. 
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